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Abstract 

In this paper we consider a probabilistic signal-to-interference-and-noise ratio (SINR) constrained 
problem for transmit beamforming design in the presence of imperfect channel state information (CSI), 
under a multiuser multiple-input single-output (MISO) downlink scenario. In particular, we deal with 
outage-based quality-of-service constraints, where the probability of each user's SINR not satisfying a 
service requirement must not fall below a given outage probability specification. The study of solution 
approaches to the probabilistic SINR constrained problem is important because CSI errors are often 
present in practical systems and they may cause substantial SINR outages if not handled properly. 
However, a major technical challenge is how to process the probabilistic SINR constraints. To tackle 
this, we propose a novel relaxation-restriction (RAR) approach, which consists of two key ingredients — 
semidefinite relaxation (SDR), and analytic tools for conservatively approximating probabilistic con- 
straints. The underlying goal is to establish approximate probabilistic SINR constrained formulations in 
the form of convex conic optimization problems, so that they can be readily implemented by available 
solvers. Using either an intuitive worst-case argument or specialized probabilistic results, we develop 
various conservative approximation schemes for processing probabilistic constraints with quadratic uncer- 
tainties. Consequently, we obtain several RAR alternatives for handling the probabilistic SINR constrained 
problem. Our techniques apply to both complex Gaussian CSI errors and i.i.d. bounded CSI errors 
with unknown distribution. Moreover, results obtained from our extensive simulations show that the 
proposed RAR methods significantly improve upon existing ones, both in terms of solution quality and 
computational complexity. 
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I. Introduction 

In multi-antenna multiuser downlinks, linear transmit beamforming has been recognized as an impor- 
tant technique, capable of leveraging quality of service (QoS) and increasing limits on the number of 
users served; see, e.g., the review article [3] and the references therein. Transmit beamforming design 
approaches developed in this context have not only been proven to offer efficient and flexible solutions 
for QoS optimization and interference management in standard downlinks, but have also been modified 
or generalized to deal with designs arising from frontier scenarios, such as relay networks [3], cognitive 
radios [4], and multicell coordinated downlinks [5]-[7]. 

In transmit beamforming, a very representative problem setting is the unicast multi-input single-output 
(MISO) downlink scenario, wherein a multi-antenna base station simultaneously transmits data streams 
to a number of single-antenna users, each stream for a designated user, by carefully directing transmit 
beams to the users. The problem of interest is to provide a signal-to-interference-and-noise ratio (SINR) 
constrained design formulation, in which transmit beamformers for the users are sought, so that each 
user is served with a QoS, characterized by the SINR, no less than a prescribed requirement, and that the 
transmit power is minimized. The SINR constrained problem is a meaningful and frequently used design 
formulation in practice, and essentially the same problem formulation can be seen in other works, such as 
those in the aforementioned frontier scenarios [6], [8], [9]. It is also a fundamentally intriguing problem. 
There are three parallel solution approaches to the problem, namely, uplink-downlink duality [10], [11], 
semidefmite relaxation (SDR) [5], [12], and the second-order cone programming (SOCP) formulation [13]. 
Each of those approaches is elegant, offering different implications both in theory and in practical 
implementations. They also serve as stepping stones for more advanced designs, such as those under 
imperfect channel state information (CSI) effects. 

The SINR constrained problem, like many other transmit optimization problems, is based on the 
assumption that the downlink CSI is perfectly available at the base station. Unfortunately, such an 
assumption generally does not hold in practice [14]. In the time division duplex (TDD) setting, where 
there is a reciprocity between the uplink and downlink channels, the downlink CSI is typically acquired 
by uplink channel estimation from training data. Channel estimation errors, which are caused by noise 
and a limited amount of training data, result in CSI errors in this setting. In the frequency division 
duplex (FDD) setting, CSI acquisition is often achieved by CSI feedback with limited rates. As a result, 
quantization errors arising from the limited feedback lead to imperfect CSI. In addition, CSI may become 
somewhat outdated if the user mobility speed happens to be faster than the CSI update speed. If one 
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uses the corrupt CSI directly to design the transmit beamformers, then the users may experience severe 
SINR outages and not be able to receive their anticipated QoS levels. 

Recently, there has been much attention on transmit beamforming designs that are robust against 
CSI errors. In particular, it is of significant interest to consider "safe" SINR constrained formulations 
under various CSI error models, where users' SINR requirements must be satisfied even with the worst 
possible CSI errors, or, alternatively, with high probability. One commonly considered formulation at 
present is the worst-case SINR constrained problem, in which the CSI errors are assumed to lie in 
a bounded set (known as the uncertainty set). This worst-case robust problem appears to be a hard 
(nonconvex) problem, since the worst-case SINR constraints are semi-infmite and indefinite quadratic. 
Several concurrent approximation schemes have been proposed to tackle the worst-case robust problem; 
notable works include the conservative SOCP formulation [15], the robust MMSE formulation [16], 
[17], and SDR [18]. The beauty of these works lies in the careful combination of robust optimization 
results [19] and problem formulations, leading to convex and tractable design solutions. 

Another safe formulation, which is the focus of this paper, is the probabilistic, or outage-based, SINR 
constrained problem. In this formulation, we assume a random CSI error model, such as the popular 
complex Gaussian model, and the SINR outage probability of each user must be kept below a given 
specification. Unfortunately, while the worst-case SINR constrained problem is considered hard to solve 
already, this is even more so with the probabilistic SINR constrained problem — Probabilistic SINR 
constraints generally have no closed form expressions and are unlikely to be easily handled in an exact 
way. Thus, one has to resort to approximate design solutions. To date, there are very few works on the 
probabilistic SINR constrained problem under the unicast downlink scenario. In [20], the authors fix the 
transmit beam directions as zero forcing and then deal with a probabilistic power control problem. In 
[21], a conservative SOCP formulation is developed using some advanced results in chance constrained 
optimization [22], [23]. A similar approach is presented in [24], where the robust MMSE formulation is 
considered. 

In this paper we propose several convex optimization solutions for approximating the probabilistic SINR 
constrained problem. Our approach is based on a relaxation-restriction (RAR) methodology. Specifically, 
in the relaxation step, we employ SDR to linearize the quadratic terms in the SINR expression. However, 
this step alone does not lead to an efficiently solvable formulation, because the probabilistic constraints 
imposed on the linearized SINR expressions are still nonconvex. We circumvent this problem in the 
restriction step, where we first derive various analytic upper bounds on the violation probability (i.e., 
the probability that the constraints on the linearized SINR expressions are violated). Such upper bounds 
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serve as sufficient conditions for the probabilistic constraints to hold, hence the term "restriction". Next, 
we show that our derived bounds are efficiently computable, which, together with the results from the 
relaxation step, leads to efficiently solvable approximations of the original probabilistic SINR constrained 
beamforming problem. It should be noted that the above restriction approach has many advantages. First, it 
allows one to generate feasible solutions to the probabilistic constraints, even when there is no closed form 
expression for the violation probability, or when the closed form expression is not efficiently computable. 
Secondly, while it may be difficult to derive closed form expressions for the violation probability, it 
is usually much easier to derive upper bounds on it, thanks to the many powerful techniques from 
the probability theory literature. Thirdly, there is usually more than one way to derive upper bounds 
on the violation probability, and this offers the possibility of trading approximation performance with 
computational complexity. These advantages will become clear in our subsequent exposition. 

The rest of this paper is organized as follows. The problem statement of the outage-based SINR 
constrained robust beamforming design problem is given in Section EL The idea of the proposed RAR 
method is introduced in Section III. In Sections IV and V, various RAR formulations for complex 
Gaussian CSI errors are developed using either robust optimization or probabilistic techniques. An RAR 
formulation for i.i.d. bounded CSI errors is also presented in Section V. Simulation results are then 
presented in Section VI, and conclusions are drawn in Section VII. 

Notations: We use boldfaced lowercase letters, e.g., a, to represent vectors and uppercase letters, 
e.g., A, to represent matrices. The notations M n , C n , S n , and HP stand for the sets of n-dimensional 
real vectors, complex vectors, real symmetric matrices and complex Hermitian matrices, respectively. The 
superscripts 'T' and k H' represent the transpose and (Hermitian) conjugate transpose, respectively. A y 
means that the matrix A is positive semidefmite. Tr(A) and A max (A) denote the trace and maximum 
eigenvalue of A, respectively. vec(A) stands for the vector obtained by stacking the column vectors of 
A. [a]i and [A]ij (or simply ctj and A%j) stand for the ith entry of a and (i,j)th entry of A, respectively. 
For a complex A, we denote by Re{ A} and Im{ A} its real and imaginary parts, respectively. I n denotes 
the n x n identity matrix. Given scalars a\, . . . , a n , we use Diag(ai, ■ ■ ■ , a n ) to denote the nxn diagonal 
matrix whose ith diagonal entry is ctj. || • || and || • represent the vector Euclidean norm and matrix 
Frobenius norm, respectively. E{-}, Prob{-}, and exp(-) denote the statistical expectation, probability 
function and exponential function, respectively. We write x ~ CJ\f(ii, C) if x — /j, is a circular symmetric 
complex Gaussian random vector with covariance matrix C >z 0. 
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II. Problem Formulation 

We focus on a downlink multiuser MISO scenario, in which the base station, or the transmitter, sends 
parallel data streams to multiple users over the same channel. The transmission is unicast; i.e., each 
data stream is exclusively for one user. The base station is equipped with Nt transmit antennae and the 
signaling strategy is transmit beamforming. Let x(t) G C * denote the multi-antenna transmit signal 
vector of the base station at time t. We have the following transmit signal model: 

K 

x(t) = ^2w k s k (t), (1) 

k=l 

where w k G C Nt is the transmit beamforming vector for user k, K is the number of users, and Sfc(i) is 
the user-/c data stream, which is assumed to have zero mean and unit power (i.e., E{|sfc(t)| 2 } = 1). It is 
also assumed that Sfe(i) is statistically independent of one another. For user i, the received signal can be 
modeled as 

yi (t) = hfx(t) + ni{t), (2) 

where hi G C^* is the channel from the base station to user i, and rii(t) is an additive noise, which is 
assumed to have zero mean and variance of > 0. 

A common assumption in transmit beamforming is that the base station has perfect knowledge of 
hi,... ,hx', i.e., the so-called perfect CSI setting. However, as discussed in detail in the Introduction, 
the base station may not have perfect CSI in general. In this work, the CSI is modeled as follows: 

hi = hi + ei, i = l,...,K, 

where hi G C 4 is the actual channel, hi G C * is the presumed channel at the base station (also 
called the imperfect CSI), and e,- L G C Nt is the respective error that is assumed to be random. Our 
development will concentrate mainly on complex Gaussian CSI errors, which is a commonly adopted 
model. Specifically, we assume that 

ei~CN{Q,Ci) 

for some known error covariance Cj y 0, i = 1, . . . , K. 

The goal here is to design beamforming vectors wi, . . . , wk such that the QoS of each user satisfies 
a prescribed set of requirements under imperfect CSI, while using the least possible amount of transmit 
power in doing so. To put this into context, let us consider users' SINRs. Under the model in (l)-(2) and 
the associated assumptions, the SINR of user i is 

SIN R, = - Iff*!' „ , 

August 5, 2011 DRAFT 



6 



To accommodate imperfect CSI knowledge at the base station, which causes uncertainties in the actual 
SINRs, we consider the following robust beamforming design problem: 

Probabilistic SINR constrained problem: Given minimum SINR requirements 71,..., jk > and 
maximum tolerable outage probabilities pi, ... , px G (0, 1], solve 

K 

min ll' u '«l| 2 (3a) 

w 1 ,...,w K ec N t z — ' 

i=l 

s.t. Prob hi „ CX{hijCi) {SINR* > 7i } > 1 - p h i = l,..., K. (3b) 

Formulation (3) is an instance of the so-called chance constrained optimization problem due to the 
presence of the probabilistic constraints (3b), and it will be the main focus of this paper. In (3), the design 
parameters p.;'s govern service fidelity, making sure that each user, say, user i, is served with an SINR 
no less than 73 at least (1 — pi) x 100% of the time. In fact, the simulation results in Section VI will 
demonstrate that a "non-robust design"; i.e., designing the beamformers by running the perfect-CSI-based 
SINR constrained problem with actual channels hi substituted by the presumed channels hi, can suffer 
from serious SINR outage. Moreover, it should be noted that there is a tradeoff between service fidelity 
and design conservatism. On one hand, it is desirable to request higher service fidelity by using small 
values with /Vs. On the other hand, the design in (3) would become more conservative as pj's decrease. 
In particular, for very small /Vs, one may end up with design solutions that have unacceptably large 
transmit power, or there may be no feasible solution to (3). 

Although the probabilistic SINR constrained problem in (3) is a meaningful design criterion, it is 
a very hard problem. The main difficulty is that the probability functions in (3b) do not yield simple 
closed form expressions for the considered CSI error distribution models. Thus, one may only resort to 
approximation methods. In the next section we will describe our proposed approximation approach. 

III. The Relaxation-Restriction Approach 

To handle the main problem (3), we propose a novel relaxation-restriction (RAR) approach. RAR-based 
methods feature the use of convex optimization techniques to approximate problem (3). Hence, they can 
be efficiently implemented by available convex optimization software. 

A. Relaxation Step 

Let us first elaborate on the first step of RAR — relaxation. The motivation is that for each i, the 
inequality SINRj > 7, is nonconvex in tci, . . . , wk', specifically, it is indefinite quadratic. We handle this 
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issue by semidefinite relaxation (SDR) [3], [25]. To illustrate SDR for the probabilistic SINR constrained 
problem, we note that problem (3) can be equivalently represented by 

K 

min VTr(Wi) (4a) 

w u ...,w K <m N * ^ 



8=1 



Pi , % = !,..., K, 

(4b) 

Wx,...,W K hO, (4c) 

rank(Wi) = 1, i = l,...,ET, (4d) 

where the connection between (3) and (4) lies in the feasible point equivalence 

Wi = Wiivf, i = l,...,K. 

The SDR of (4) works by removing the nonconvex rank-one constraints on Wi; i.e., to consider the 
relaxed problem 

K 

min VTr(Wi) (5a) 

Wi,...,w K ea N i 



s.t. Prob < 



(hi + ei) H f -Wi -TwA (hi + e-) > a? \ > 1 



Pi, i = l,...,K, 



(5b) 

W U ...,W K h 0. (5c) 

The merit of this relaxation is that the inequalities inside the probability functions in (5b) are linear in 
W\, . . . , Wk, which makes the probabilistic constraints in (5b) more manageable. An issue that comes 
with SDR is the solution rank — the removal of rank(Wj) = 1 means that the solution (Wi, . . . , Wk) 
to problem (5) may have rank higher than one. We shall come back to this issue after presenting the 
restriction step. 

B. An Information Theoretic Interpretation of the Relaxation Step 

The SDR problem (5) has an alternative interpretation from an information-theoretic point of view. 
Here, we briefly describe this interpretation and the resulting implications before moving to the restriction 
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step of RAR. Consider a general transmission model: 

K 

x(t) = ^2xi(t), 

i=l 

where Xi(t) £ C Nt is the transmit signal intended for user i. Note that in contrast to the original transmis- 
sion model in (1), where we fix the transmit scheme as beamforming by setting Xi(t) = WiSi(t), here we 
do not assume any specific transmit structure. As a slight abuse of notations, let Wi = K{xi(t)xf (t)} be 
the transmit covariance corresponding to user i. From an information-theoretic perspective, the achievable 
rate of user i may be formulated as 

hfWjhj 

where the rates in (6), in bits per channel use, are achieved when Xi(t) are Gaussian distributed (i.e., 
Gaussian codebook). One can easily verify that the SDR problem (5) is precisely the following rate 
optimization problem: 

K 

min VTr(Wi) (7a) 
Wi,...,w K ea N t ^ 
i=i 

s.t. Prob^cj^^.) {Ri > log 2 (l + 7i)} > 1 - p u i = 1,. . . ,K, (7b) 
W 1 ,...,W K hO. (7c) 

Specifically, the SDR equivalent (7) is a total power minimization problem that aims to ensure that each 
user is served with a minimal rate of log 2 (l + 7«) bits per channel use, with an outage probability no 
greater than pi. 

With this interpretation of the SDR, we can deduce an interesting implication: If the SDR solution 
(Wi, . . . , Wk) to (5) does not yield a rank-one structure, then an alternative to transmit beamforming is 
to find another practical physical-layer scheme — e.g., a space-time code with appropriate precoding — 
to adapt the transmit structures stipulated by the transmit covariances (W\, . . . , Wk)- While our main 
interest in this paper is still in transmit beamforming, it is worthwhile to keep such a parallel possibility 
in mind, since it eliminates the need for rank-one transmit covariances. 

C. Restriction Step 

Let us continue to illustrate the second step of RAR — restriction. The relaxation step alone does not 
provide a convex approximation of the main problem (3). The SDR probabilistic constraints (5b) remain 
intractable, although they appear to be relatively easier to handle than the original counterparts in (3b). 
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The restriction step aims to find a convex approximation of (5b), in a restrictive or conservative sense. 
More precisely, in the context of the probabilistic SINR constrained problem, the restriction step entails 
finding a solution to the following: 



Challenge 1: Consider the following chance constraint: 

Pmb{e H Qe + 2Re{e H r} + s > 0} > 1 - p, (8) 

where e G C n is a standard complex Gaussian vector (i.e., e ~ CM(0, I n )), the 3-tuple (Q,r,s) G 
HP x C n x M is a set of (deterministic) optimization variables, and p G (0, 1] is fixed. Find an efficiently 
computable convex restriction of (8); i.e., find an efficiently computable convex set S C HP xPxixC' 
such that whenever (Q, r, s, t) G S, the 3-tuple (Q, r, s) G HP xPxl is feasible for (8). 



Note that in the construction of the convex set S, we are allowed to include an extra optimization 
variable t G C £ , in addition to the original optimization variables (Q, r, s) G HP x C x R. Although the 
precise role of t will depend on how the set S is formulated, it suffices to think of t as a slack variable. 

It is not hard to see that the SDR probabilistic constraints in (5b) fall in the scope of Challenge 1. 
Indeed, for each constraint in (5b), the following correspondence to (8) can be shown: 

Q = C}' 2 (±Wi - £ W k ) C)l\ r = C]' 2 (im - £ W k ) hi, (9a) 

\ 7? / \ 7j fc* / 

s = hf I - Wi - ^ W k ] ^ - ^, p = ^. (9b) 

V 7 ' / 

The development of convex restriction methods plays a crucial role in RAR, and this will be our 
focus in subsequent sections. By replacing each probabilistic constraint in (5b) with a convex restriction, 
we will obtain a convex approximation of the original probabilistic SINR constrained problem. Table I 
summarizes all the RAR methods to be proposed in later sections. Note that each RAR method is based 
on a different convex restriction. Moreover, all the RAR formulations in Table I are conic problems with 
linear matrix inequality constraints and/or second-order cone constraints, which can be easily solved by 
off-the-shelf convex optimization software [26]. 

The last step of RAR is to provide a feasible beamforming solution (w\, . . . , wk) to the main 
problem (3) by using the RAR solution (Wi, . . . , Wk)- As is common in all SDR-based methods, 
the Wj's obtained from RAR may have rank higher than one. A standard way of tackling this issue is 
to apply some rank-one approximation procedure to (Wi, . . . , Wk) to generate a feasible beamforming 
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TABLE I. The proposed RAR formulations. 



Method 


RAR Formulation 


Method I: 
Sphere bounding 

(for complex Gaussian 
CSI errors) 


K 

min VTr(Wi) 

Wj eu N t em, ' 

i=l,...,K 1=1 

Q t + Ul Nt r, ^ (10) 
s.t. h 2 bO, i = l,...,K, 

Wi,...,W K hO, ti,...,t K >0; 

where Qi, n and s, are defined in the same way as (9), and dt — ./$ _ 2 1 (1 — Pi)/2, i 
1,...,K. 


Method II: 

Bernstein-type 

inequality 

(for complex Gaussian 
CSI errors) 


min VTr(Wi) 

i=l,...,K l ~ 

s.t. Tr(Qi) - y/-2 ln(pi) • x, + ln(pi) ■ y< + s, > 0, i = 1, . . . , K, 

fvec(Q0l ^ . 1 (ID 

2/iIjv t + Qt h 0, i = 1, . . . , K, 
Vi,..., VK > 0, W U ...,W K h 0; 
where Qi, n and s* are defined in the same way as (9), i — 1, . . . , K. 


Method III: 
Decomposition into 
independent parts 

(for complex Gaussian 
CSI errors) 


K 

min VTr(Wi) 
i=i,...,K 1_ 

s.t. + Tr(Qi) > 2^- ln( Pl ) • (a* + * = 1. ■ • • » K 
■^\\ri\\ < Xi, i = l,...,K, 
Vi [|vec(Qi)]| < i/i, i = 1, 
>: 0; 

where Qi, ri and Sj are defined in the same way as (9), Qi is chosen such that 6i + ln(l — Si) = 
ln(pi), and Vi = \J- ln(pi)/0f, i = 1, . . . ,K. 


Method IV: 
Decomposition into 
independent parts 

(for elementwise i.i.d. 
and bounded CSI errors 
with mean and vari- 
ance <7e, but otherwise 
unknown distribution) 


K 

min VTr(Wi) 

W; EM T ,i~ EJK f r ... 

i=l,.!.,K 1=1 

s.t. si + a 2 e ■ Tr(Qi) > 2V-ln(p») ■ ESiM«> Vi > 

V2||ri|| < [U}o, i = l,...,K, (13) 

£ = l,...,2N t , i = l,...,K, 
Wx,...,W K t 0; 

where Qi, r» and Si are defined in the same way as (32), At are defined in the same way as 
Table II, and Vjj = l/y/8 and Vjk = 1 if j / fc. 
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solution . . . ,wk) to (3); see [25] for a review and references. In our setting, we apply a Gaussian 
randomization procedure to a non-rank-one RAR solution. The procedure is provided in Algorithm 1; 
the spirit follows that of [27], and readers are referred to [27] for an exposition of the idea. We should 
point out that obtaining a feasible RAR solution does not imply that we can always generate a feasible 
solution to the main problem (3). This issue has also been identified before in the context of multigroup 
multicast beamforming with perfect CSI [27]. However, if the RAR solution happens to give rank- 
one Wi for all i, then we can simply solve the rank-one decomposition Wi = WiW^ and output the 
corresponding (wi, . . . , wr) as the approximate beamforming solution. For such instances, it can be 
easily verified that (w\, . . . ,wk) is already feasible for the main problem (3). Rather surprisingly, we 
found that the proposed RAR methods returned rank-one solutions in almost all the simulation trials we 
ran 1 . Such empirical finding provides another interesting implication when we consider the information- 
theoretic interpretation in the last subsection: Since the RAR methods are essentially the same as convex 
restrictions of the outage-based rate optimization problem in (7), the numerical observation that RAR 
solutions are almost always rank-one somehow hints that transmit beamforming may inherently be an 
optimal physical-layer scheme, at least for the outage-based unicast multiuser MISO downlink scenario 
considered here. 

As a summary to the solution approximation aspect discussed above, in most cases a simple rank-one 
decomposition Wi = Wiwf 1 of the RAR solution suffices to produce a feasible solution (wi, . . . , wk) to 
the main problem (3). The more complicated solution approximation procedure in Algorithm 1, proposed 
for instances where the RAR solution is not of rank one, is rarely needed in our empirical experience. 

IV. RAR Method I: Sphere Bounding 

In this section we describe our first convex restriction method for Challenge 1. The method is based 
on two key ingredients. The first is the following lemma: 

Lemma 1 Consider Challenge 1. Suppose that we have a set B C C n that satisfies 

Prob{e e B} > 1 - p. (14) 

Then, the following implication holds: 

S H Q8 + 2Re{5^r} + s > 0, 

=>■ Eq. (8) in Challenge 1 holds. (15) 

for all 8 e B 

'A similar phenomenon was observed in a different problem setting, namely that of the worst-case SINR constrained 
design [18]. 
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Algorithm 1: Gaussian randomization procedure for RAR 
Given : A number of randomizations L, and an optimal solution (W*, . . . , W£) to an employed 



1 for £ = 1, ... ,L do 

2 generate random vectors ~ CA/"(0, W*), i = 1, . . . , K ; 

3 set beam directions uf' = w\ ' /\\w\' \\, i = 1, . . . , K; 

(£) (£) (£) (£) 

4 let p\ , . . . ,p K be beam powers and obtain p\ , . . . ,p K as follows: substitute 

W~i = piVLp {v!p) H , i = 1, . . . , K, into the RAR problem, solve the problem with respect to 
Pi, ■ ■ ■ ,pk > 0, and set p(\ . . . as its solution if the problem is feasible; also, set to 
be the associated optimal objective value if the problem is feasible; otherwise set = oo ; 

5 end 

6 £* = arg min^ = i j ^ 

Output: w* = y pf ^uf \ i = 1, . . . , K, as an approximate solution to the main problem (3). 



The proof of Lemma 1 is simple and is given as follows. Let p(e) denote the probability density function 
of e. Suppose that (14) and the left-hand side (LHS) of (15) hold. Then, we have the following chain: 



Hence, Eq. (8) is satisfied. 

Lemma 1 suggests that we can approximate the chance constraint in (8) in a conservative (or restrictive) 
fashion by using the worst-case deterministic constraint on the LHS of (15). Moreover, it can be easily 
seen that the same idea applies to general chance constraints; i.e., the quadratic functions in (8) and (15) 
may be replaced by any arbitrary function. Such an insight (i.e., using worst-case deterministic constraints 
to approximate (general) chance constraints) have been alluded to or used in many different contexts; 
e.g., [28], [29] in optimization. Here, we are interested in the chance constraint in (8), which involves 
a quadratic function of the standard complex Gaussian vector e. In our method, we choose 6 to be a 
spherical set; i.e., 



RAR formulation. 




B = {SeC n | ||«|| <d} 



(16) 
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where d is the sphere radius. It can be shown that by choosing 

where $>~2 (•) is the inverse cumulative distribution function of the (central) Chi-square random variable 
with m degrees of freedom, Eq. (14) is satisfied. 

The second ingredient is the so-called 5-lemma, which enables us to turn the infinitely many constraints 
on the LHS of (15) into a set of tractable constraints. The 5-lemma is given as follows: 

Lemma 2 ( 5-lemma [30] ) Let fi(x) = x H QiX + 2Re{x H Ti} + Sj for i = 0,1, where x £ C n and 
(Qi,Ti, Si) E I" x C n x M.for i = 0, 1. Suppose that there exists an x G C n satisfying f\(x) < 0. Then, 
the following statements are equivalent: 

1- fo(x) > Ofor all x € C n satisfying f x (x) < 0. 

2. There exists a t > such that 



Qo 


ro 


+ t 


Qi 


ri 


. 




/ 1 




so 




si 



By the 5-lemma, the LHS of (15), with B given by (16), can be equivalently represented by an LMI of 
the form (17), where (Qo, To, «o) = {Qi r i s ) an d (Qi, fx, «i) = (In: 0, —d 2 ). We therefore have built a 
convex restriction for Challenge 1. To summarize, we have the following: 



Method I for Challenge 1 (Sphere bounding): The following feasibility problem is a convex restriction 


of (8) in Challenge 1: 






Find Q,r,s,t 






Q + tl n r 


to, 


s.t. 




r H s - td 2 


t > 0, 




where d= J&~} (l-p)/2. 

V X2ti 







By first applying SDR and then Method I to the probabilistic SINR constrained problem (3), we obtain 
the RAR formulation (10) in Table I. Interestingly, this formulation turns out to be similar to that of 
the worst-case robust SDR problem considered in [18]. However, it should be noted that the prior work 
does not consider outage probability constraints. Moreover, we show a way of using the worst-case 
robust formulation to deal with the probabilistic SINR constrained problem. Finally, by incorporating the 
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bisection scheme proposed in [21], which will be considered in our simulations in Section VI (Example 
3), we will be able to further improve the performance of the sphere bounding RAR method. 

V. Probability Inequality Approaches 

The reader may notice that the development of Method I is strongly motivated by the worst-case robust 
optimization paradigm. Indeed, the problem on the LHS of the implication (15) is precisely a robust 
feasibility problem with uncertainty set B. By choosing B judiciously, it is shown that the violation 
probability Probje^Qe + 2Ke{r H e} + s < 0} can be controlled, and the resulting robust feasibility 
problem is a convex restriction of (8). However, this approach has an intrinsic drawback, namely, it is 
difficult to define and analyze an uncertainty set B other than those that have very simple geometry, such 
as the spherical set considered in the previous section. Consequently, it is not clear whether there exist 
other choices of B that would lead to better convex restrictive approximations. 

As it turns out, one can circumvent the above drawback by using analytic upper bounds on the violation 
probability to construct efficiently computable convex restrictions of (8). Specifically, suppose that we 
have an efficiently computable convex function f(Q,r,s,t), where t is an extra optimization variable, 
such that 

Prob{e^Qe + 2Re{e H r} + s < 0} < f(Q, r, s, t). (18) 

Then, the constraint 

f(Q,r,s,t)<p (19) 

is, by construction, a convex restriction of (8). An upshot of this approach is that there are many available 
techniques for constructing such upper bounds, and each of those bounds yields a convex restriction of (8). 
Moreover, it is known [19, Chapter 4] that under some fairly mild conditions, every convex restriction 
corresponds to a robust feasibility problem with a suitably defined uncertainty set. Thus, the above 
approach can be viewed as an enhancement of Method I, in the sense that it provides a handle on more 
sophisticated uncertainty sets that are difficult to construct directly. 

A. Method II: Bernstein-Type Inequality 

Let us now illustrate the above approach by showing how a Bernstein-type inequality for Gaussian 
quadratic forms can be used to construct a convex restriction of (8). Our approach relies on the following 
lemma due to Bechar [31]: 
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Lemma 3 Let e ~ CJ\f(0, I n ), Q e BP and r £ C n . 77ie«, /or any rj > 0, we have 

Prob{e H Qe + 2Re{e^r} > r(r/)} > 1 - e"* 7 , (20) 
where the function T : M++ — > M. is defined by 

T{ v ) = Tr(Q) - IIQIH + 2||r P - V X+(Q), (21) 

with X + (Q) = max{A max (-Q),0}. 

Lemma 3 is obtained by extending the corresponding result in [31] for quadratic forms of real-valued 
Gaussian random variables. The inequality in (20) is a so-called Bernstein-type inequality 2 , which bounds 
the probability that the quadratic form e H Qe+2Ke{e H r} of complex Gaussian random variables deviates 
from its mean Tr(Q). 

Since T(rj) is monotonically decreasing, its inverse mapping T _1 : M — > M++ is well defined. In 
particular, the Bernstein-type inequality in (20) can be expressed as 

Prob{e^Qe + 2Re{e H r} + s > 0} > 1 - e" T_1( " s) . (22) 

As discussed in (18) and (19), the constraint e~ T < p, or equivalently, 

- V~ 2 Hp)\/ WQWl + 2 IMI 2 + Hp) ■ + s > o (23) 

serves as a sufficient condition for achieving (8). 

While it is not obvious at this stage whether (23) is convex in (Q, r, s) or not, a crucial observation 
is that (23) can be equivalently represented by the following system of convex conic inequalities: 

Tr (Q) - V"21n(p) • h + ln(p) • t 2 + s > 0, (24a) 

^IIQIlF + 2||r|| 2 < h, (24b) 

t 2 In + QhO, (24c) 

t 2 > 0, (24d) 

where t\,t 2 £ M are slack variables. Therefore, formulation (24) is an efficiently computable convex 
restriction of (8). We now summarize the Bernstein-type inequality method as follows: 

2 Roughly speaking, a Bernstein-type inequality is one which bounds the probability that a sum of random variables deviates 
from its mean. The famous Markov inequality, Chebyshev inequality and Chernoff bounds can all be viewed as instances of 
Bernstein-type inequalities. 
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Method II for Challenge 1 (Bernstein-type inequality method): The following feasibility problem is 
a convex restriction of (8) in Challenge 1: 

Find Q,r,s,t 
s.t. Tr (Q) - V-21n(p) • h + ln(p) • t 2 + s > 0, 
^/||Q||| + 2||rP <t 1? 
t 2 In + QhO, 
t 2 > 0. 

Upon applying Method II to (5), we obtain the RAR formulation (11) in Table I. As can be easily 
seen from the formulations (10) and (11), the latter has a more complex constraint set and thus a higher 
computational complexity in general. However, it will be shown later that the Bernstein-type inequality 
method (11) exhibits better approximation performance than the sphere bounding method. 

B. Method III: Decomposition into Independent Parts 

For both the sphere bounding and Bernstein-type inequality methods, the resulting convex restrictions 
of (8) contain linear matrix inequality constraints. As such, they could be computationally costly when the 
problem size is large. It turns out that one can also develop a convex restriction of (8) that contains only 
second-order cone constraints. The resulting formulation can thus be solved more efficiently than those 
developed using the sphere bounding or Bernstein-type inequality method. The idea is to first decompose 
the sum e H Qe + 2Re{e^r} + s into several parts, each of which is a sum of independent random 
variables. Then, one bounds the moment generating function of each of those parts and stitch the results 
together to obtain an analytic upper bound on the violation probability [32]. To illustrate this approach, 
let Q = UAU H be the spectral decomposition of Q, where A = Diag(Ai, . . . , A n ) and Ai, . . . , A n are 
the eigenvalues of Q. Since e ~ CAf(0, I n ) and U H is unitary, we have e = U H e ~ CAf(0, I n ). Thus, 
we can write 

* = e H Qe + 2Re{e H r} = e H Ae + 2Re{e H r} = ^ g + % 
Now, observe that both 

n n 

^ q = e H Ae = ^2\ e \e e \ 2 and ^ = 2Re{e H r} = 2 ^ (Re{r^}Re{e^} + Im{r*}Im{e*}) 



August 5, 2011 



DRAFT 



17 



are sums of independent random variables. Moreover, it can be shown that for any fixed 6 < 1, 

E {exp (9(\e e \ 2 - 1)) } = < exp (v 2 9 2 ) , 

E{exp(0- 2Re{e*})} = E{exp((9 • 2Im{e e })} = exp ( -6 2 J for 6» G M, 



s 2 

where u = (— (0 + In (l — 0)) /6 2 ) < oo. Thus, for any p±,p2 > such that pi + P2 = 1, the chain 
of inequalities 

E{exp(u(Tr(A) - *))} = E jexp ( Pl ■ ^(^g " Tr(A)) + p 2 ■ } 

< Pl E|exp^-^(^-Tr(A))^| +p 2 E|exp^-^^| (25) 



Pi 



=i 



Pi 



~[E\exp(--\ e (\e e \ 2 -l))\ (26) 



< 



+ p 2 f[ E jexp (-^-2Re{r e }Re{e e }^ j E jexp ^-^2Im{r^}Im{e^}^ | 



pi exp 



(27) 

holds whenever —u\( jp\ < 6 for £ = 1, . . . , n, where (25) follows from Jensen's inequality and (26) 
follows from the independence of the random variables in ^ q and By setting 

t=l 

we see from (27) that the inequality 

E{exp(u(Tr(A) - *))} < pi • exp (u 2 T 2 ) +p 2 ■ exp (u 2 T 2 ) = exp (u 2 T 2 ) 
holds whenever \u\T < 6v. In particular, by Markov's inequality, it can be shown that for any ( > 0, 
Prob{Tr(A) - ^ > (} < inf ( exp(-uC) • E {exp(u(Tr(A) - *))} ) 

0<u<9v/T l J 



exp ( + (0t>) 2 ) forC>2^T. 



exp ( 1 



(28) 



I 1 

Now, set C = s + Tr(A). Then, the LHS of (28) becomes 

Prob{^ + s < 0} = Prob {e H Qe + 2Re{e^r} + s < 0} . 
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In particular, by imposing the constraint that the right-hand side of (28) is less than p and using the fact 
that Tr(A) = Tr(Q) and J21=i = IIQIIf> followed by some tedious derivations (see [32] for details), 
we obtain the following method for Challenge 1. 



Method III for Challenge 1 (Decomposition into Independent Parts): Given a parameter 9 < 1, let 

/ 9 + \n(l-e)\ 

v = — 



and 



e 2 



2^Hp), if 9v > v /Z H^), 



9v otherwise. 



Hp) 

Ov ■ 

Then, the following feasibility problem is a convex restriction of (8) in Challenge 1: 

Find Q,r,s,t 
s.t. s + Tr (Q) - vih + h) > 0, 

v\\Q\\f < h- 



Observe that in Method III, we have the flexibility to choose the parameter 9. Ideally, 9 should be chosen 
so that both p and v are small, since then the constraints in (29) are easier to satisfy. However, as can be 
seen from the definition, /i and v cannot be chosen independently of each other. Our simulation results 
suggest that it is better to have a smaller value of /x; i.e., choose 9 (and hence v) such that fi = 2^J — \n(p). 
Specifically, for a given p G (0, 1), we choose 9 such that v is minimized and p = 2^/— ln^). This can 



be achieved by solving 9v = \J — fn(p), or equivalently, 

9 + ln(l - 9) = ln(p), (30) 

which can be done numerically. We remark that for small values of p (say, p 6 (0, 0.2)), the solution 9 
to (30) can be approximated by 

0^1- exp(ln(p) - 1). 

C. Variation on a Theme: i.i.d. Bounded CSI Errors with Unknown Distribution via the Decomposition 
Approach 

An advantage of the decomposition approach outlined above is that it can be applied to cases where 
the distribution of the random vector e is not Gaussian. As an illustration, let us generalize the setting 



August 5, 2011 



DRAFT 



19 



considered in the previous section and develop an RAR method for handling the elementwise i.i.d. bounded 
support model with unknown distribution. In this model, the real and imaginary parts of the CSI error 
vector Bi are assumed to be independent and have i.i.d. components. Each component has zero mean and 
is supported on, say, [— e;, ej, where > 0. Again, we pose the restriction step in RAR as the following 
generic challenge: 



Challenge 2: Consider the following chance constraint: 

Prob{e T Qe + 2e T r + s > 0} > 1 - p, (31) 

where e E W n is a mean-zero random vector supported on [— \/3, V3] n with independent components, 
the 3-tuple (Q, r, s) G S n x W 1 x 1 is a set of optimization variables, and p G (0, 1] is fixed. Find 
an efficiently computable convex restriction of (31). 



Q 



Note that Challenge 2 and the SDR probabilistic SINR constrained problem (5) are related via the 
following identification: 

Re^Wi-Yk^Wk) hi}' 
Juifywi-Ek&Wk) hi}_ 



V3 



(32a) 



(32b) 



hf \-Wi-Y,W k \ hi-ol 



P = Pi 



(32c) 



To tackle Challenge 2 using the decomposition approach, we first observe that the sum * = e T Qe+2e T r 
can be written as 

n n 
i=l l<£^j<n 1=1 



1=1 

II 



•) ) + ( 2 Qik e i e k 



+ 2 ^ 



a e 2 .TV(Q) + J]^ + ^, 



where o\ = E{ef} and the sets Ai, ■ ■ ■ ,A n are defined as in Table II. In other words, if the (J, k)th 
entry of the table is labeled At, then (j, k) G Ae- 

Using Table II, it is not hard to verify that each of the terms ^ q ±, ^ q 2, ■ ■ ■ , ^qn, is a sum of 
independent random variables. Thus, by bounding their moment generating functions and using an 
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TABLE II. Construction of the sets A\, . . . ,A> 





1 


2 




n - 1 


n 


1 


Ai 


A 2 




An-l 


An 


2 


A2 


As 




An 


Ai 














n- 1 


An-l 


An 




A„-3 


An-2 


n 




Ai 




A n -2 


An-l 



argument similar to that in the previous subsection, we obtain the following method for Challenge 2 
(again, see [32] for details): 



Method IV for Challenge 2 (Decomposition into Independent Parts): 


The following feasibility 


problem is a convex restriction of (31) in Challenge 2: 




Find Q,r,s,t 




n 

s.t. s + al ■ Tr(Q) > 2sJ - ln(p) • ^ t e , 

e=o 




V2\\r\\ < t , 




( V /2 




v )kQ]k <U, £ = l,...,n, 




\mzAt j 




where ai = E{ef }, Vjj = 1/V& and vjk = 1 if j 7^ k, for (j, k) S At and I 


= 1,. . . ,n. 



VI. Simulation Results 

This section shows an extensive set of simulation results illustrating the performance of the proposed 
RAR methods. 

Let us first describe the general simulation settings. We employ a universal QoS specification for all 
users; i.e., 71 = • • • = jk = 7> Pi = • • • = Pk — P- The users' noise powers are identical and fixed 
at a\ = • • • = a\ = 0.1. In each simulation trial, the presumed channels {hi]f =1 are randomly and 
independently generated according to the standard complex Gaussian distribution. 

Next, we provide some implementation details of the RAR methods. The RAR problems (those in Table 
I) are solved by the conic optimization solver SeDuMi [33], implemented through the now popularized 
and very convenient parser software CVX [26]. Then, we check whether a solution (W\, . . . , Wk) to an 
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RAR problem is of rank one or not. If yes, then the rank-one decomposition, Wi = Wiwf 1 Vi, is used 
to obtain a beamforming solution (w\, . . . ,wk)- Otherwise, the Gaussian randomization procedure in 
Algorithm 1 is called to generate a feasible (wi, . . . , w^)- Numerically, we declare that (W\, . . . , Wk) 
is of rank one if the following conditions hold: 

A ;7^ } > 0.99 forain = l,...,K; 

i.e., the largest eigenvalue of Wj is at least 100 times larger than any of the other eigenvalues. Moreover, 
we say that an RAR method is feasible if the RAR problem has a feasible solution and the subsequent 
beamforming solution generation procedure is able to output a feasible (w\, . . . , wk)- 

The RAR methods are benchmarked against the probabilistic SOCP methods in [21]. The latter are 
also implemented by SeDuMi through CVX. To provide a reference, we also run a conventional perfect- 
CSTbased SINR constrained design (e.g., [13]), where the presumed channels {hi}f =1 are used as if 
they were perfect CSI. We will call this the "non-robust method", for convenience. 



A. Simulation Example 1 

We start with the simple case of N t = K = 3; i.e., three antennae at the base station, and three users. 
The CSI errors are spatially i.i.d. and have standard complex Gaussian distributions; i.e., C\ = ■ ■ ■ = 
Cr = a elN t > where a\ > denotes the error variance. We set o\ = 0.002. The SINR requirement is 
7 = lldB. The outage probability requirement is set to p = 0.1, which is equivalent to having a 90% or 
higher chance of satisfying the SINR requirements. 

First, we are interested in examining the actual SINR satisfaction probability, ProbjSINRj > 7}, of 
the various methods. Figure 1 shows the histograms of the actual SINR satisfaction probabilities over 
different channel realizations. To obtain the histograms, we generated 500 realizations of the presumed 
channels {hi}f =1 . Then, for each channel realization, the actual SINR satisfaction probabilities of all 
methods were numerically evaluated using 10, 000 randomly generated realizations of the CSI errors 
{ei}f =l , which should be sufficient in terms of the probability evaluation accuracy. Figure 1 validates 
that our RAR methods (and the existing probabilistic SOCP method) indeed adhere to the 90% SINR 
satisfaction specification. There are two interesting observations, as can be seen from the figure. The first 
is with the non-robust method. While the non-robust method is, by nature, expected to violate the SINR 
outage specification, its actual SINR satisfaction probabilities are below 50% for most of the channel 
realizations, which is severe. This reveals that the perfect-CSTbased design can be quite sensitive to 
CSI errors. The second is with the conservatism of the various robust methods. The probabilistic SOCP 
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channel realizations 



| Probabilistic SOCP 



Counts over different 100 _ 
channel realizations 
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channel realizations 
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channel realizations 
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Actual SINR satisfaction probability (Mean = 0.38043) 



Actual SINR satisfaction probability (Mean = 1) 



Actual SINR satisfaction probability (Mean = 0.99526) 



Actual SINR satisfaction probability (Mean = 0.96695) 



Actual SINR satisfaction probability (Mean = 0.99559) 
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0.9 1 




0.9 1 



Fig. 1. Histograms of the actual SINR satisfaction probabilities of the various methods. N t = K = 3; C\ 
C K = 0.002Jjv t ; 7 = 11 dB; p = 0.1. 



method has its actual SINR satisfaction probabilities concentrating at 100%, which indicates that it may 
be playing too safe in meeting the outage specification. By contrast, our RAR methods seem to be less 
conservative. Particularly, among the three methods, RAR Method II (Bernstein-type inequality) appears 
to be the most relaxed as observed from its histogram. 

Next, we investigate the conservatism of the various robust methods by evaluating their feasibility 
rates; i.e., the chance of getting a feasible beamforming solution under different channel realizations. 
Similar to the last investigation, 500 channel realizations were used. The obtained result is shown in 
Figure 2(a), where the feasibility rates of the various methods are plotted against the SINR requirements 
7. Remarkably, the three RAR methods yield feasibility rates much higher than that of the probabilistic 
SOCP method. In particular, RAR Method II has the best feasibility rate performance, which is consistent 
with the SINR satisfaction probability result we noted in Figure 1 . The feasibility rates of RAR Methods 
I and III are a close match: For 7 > 9dB, RAR Method I slightly outperforms RAR Method III; for 
7 < 9dB, we see the converse. 

In addition to the feasibility rate, it is important to examine the transmit power consumptions of the 
design solutions offered by the various robust methods. Figure 2(b) shows the result. It was obtained based 
on channel realizations for which all methods yield feasible solutions at 7 = lldB; 181 such realizations 
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were found out of 500 realizations (the same realizations used in the last result in Figure 2(a)). As can 
be seen from Figure 2(b), RAR Method II yields the best average transmit power performance, followed 
by RAR Methods I and III (with Method I exhibiting noticeably better performance for 7 > 15dB), and 
then the probabilistic SOCP method in [21]. As a reference, we also plot the transmit powers of the 
non-robust method in the figure, so as to get an idea of how much additional transmit power would be 
needed for the robust methods to accommodate the outage specification. We see that for 7 < lldB, the 
transmit power difference between an RAR method and the non-robust method is about 1.5dB, which 
is reasonable especially when compared to the probabilistic SOCP method. The gaps gradually widen, 
otherwise. This seems to indicate that imperfect CSI effects are more difficult to cope with when we 
demand higher SINRs. 

Now, let us consider the computation times of the various robust methods. The result is illustrated 
in Figure 3. To obtain this result, we use a desktop PC with 2.13GHz CPU and 3GB RAM. Moreover, 
instead of calling the convenient parser CVX, we use direct SeDuMi implementations of all the methods, 
done by careful manual problem transformation and programming. The reason of doing so is to bypass 
parsing overheads, which may result in unfair runtime comparisons. From the figure, we see that the 
runtime ranking, from the shortest to longest, is: RAR Method III, RAR Method I, RAR Method II, and 
the probabilistic SOCP method. Interestingly and coincidently, the runtime ranking of the RAR methods 
is exactly the opposite of their performance ranking we see in the previous simulation result. 

As the last result in this example, we numerically inspect a technical issue that has much implication to 
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N t = K 



Fig. 3. Average runtimes of the various methods. 



TABLE III. Ratios of rank-one RAR solutions. 



p 


0.1 


0.01 


7 (dB) 


3 


7 


11 


15 


3 


7 


11 
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Method I 


464/464 


448/448 


404/404 


292/292 


450/450 


424/424 


343/343 


225/225 


Method II 


489/489 


475/475 


441/441 


363/363 


479/480 


463/463 


428/428 


322/322 


Method III 


488/488 


453/453 


389/389 


267/267 


476/476 


421/421 


306/306 


144/144 



the RAR approach — how frequent do the RAR problems yield rank-one solutions. Recall that rank-one 
RAR solution instances have the benefits that the beamforming solution generation is simple (simple 
rank-one decomposition, no Gaussian randomization), and that feasibility of the RAR problem directly 
implies that of beamforming solution generation. Table III shows the result. There is a ratio in each field. 
The denominator is the realizations count for which the RAR problem is feasible, and the numerator 
is the realizations count for which the RAR problem yields a rank-one solution. Again, 500 channel 
realizations were used. Curiously, almost all the fields in Table III indicate rank-one solution all the time. 
We encountered only one non-rank-one instance out of 480 for the setting of p = 0.01, 7 = 3dB, RAR 
Method II. We therefore conclude, on the basis of numerical evidence, that occurrence of high-rank RAR 
solutions is very rare for the unicast outage-based SINR constrained problem considered here. 
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Fig. 4. Performance under spatially correlated Gaussian CSI errors. N t = K = 8; p = 0.01; er^ = 0.002. 
B. Simulation Example 2 

This example considers more challenging settings, described as follows: N t = K = 8; spatially 
correlated CSI errors where C\ = ■ ■ ■ = Ck = C e , 



and <Tg = 0.002; p = 0.01 (or 99% SINR satisfaction probability). We do not run the probabilistic SOCP 
method in [21], since, as seen in Figure 3, it is computationally very demanding for large problem sizes. 
The same simulation method in Simulation Example 1 was used to produce the results here. Figure 4 
shows the resulting feasible rates and average transmit powers. A minor simulation aspect with the transmit 
power performance plot in Figure 4(b) is that we choose 7 = 7dB as the pick-up point of feasible channel 
realizations of all the methods. We can see that, once again, RAR Method II offers superior performance 
over the others. Another observation is that RAR Method III manages to outperform RAR Method I this 
time. 

Figure 5 illustrates another set of results, where we increase the CSI error variance a\ from 0.002 to 
0.01. The number of users is set to K = 6. The feasible realizations pick-up point is 7 = 13dB. We can 
see similar performance trends as in the previous result in Figure 4. 

C. Simulation Example 3 

One might have noticed from Simulation Example 1, Figure 1, that none of the robust methods 
yield actual SINR satisfaction probabilities at 1 — p. This means that the robust methods are, to certain 




m,n — 
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Fig. 5. Performance under spatially correlated Gaussian CSI errors. Nt = 8; K = 6; p = 0.01; a 2 e = 0.01. 



extent, conservative. If more computations are allowed, this conservatism may be mitigated by running 
a bisection scheme. Such a scheme was first proposed in [22] in the context of chance constrained 
optimization, and then was adopted in [21] for the probabilistic SOCP method. The idea is to fine tune 
some design parameters relevant to the outage requirement. If a design solution is found to satisfy the 
outage specification well, then we adjust the design parameters to relax the outage requirement (e.g., for 
RAR Method I, decreasing d\, . . . , Ak) and rerun the design problem. Otherwise, we do the opposite 
and rerun the design problem. The above step is done repeatedly following a bisection search, requiring 
the design problem to be solved multiple times. The bisection search also requires a validation procedure 
for satisfiability of the outage specification, which can be done using the Monte-Carlo based validation 
procedure in [22]. It is clear from the above discussion that the bisection scheme can also be applied 
to all the RAR methods. For more complete descriptions of the bisection scheme in the context of the 
probabilistic SINR constrained beamforming problem, readers are referred to [1], [2], [21]. 

Figure 6 shows how bisection may improve the performance. The simulation settings are N t = K = 5, 
p = 0.1, C\ = ■ ■ ■ = Ck = 0.002i"7v t > and the feasible realizations pick-up point at 7 = 9dB (for the 
transmit power performance evaluations only). We can see that all the robust methods, after applying 
the bisection scheme, exhibit improved performance. Notwithstanding, we also see that RAR Method II 
without bisection already gives performance quite on a par with the bisection-aided methods. 
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Fig. 6. Performance with bisection. Nt = K = 5; p = 0.1; spatially i.i.d. Gaussian CSI errors with a"l = 0.002. 

TABLE IV. Ratio of rank-one RAR solutions. 
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D. Simulation Example 4 

This example demonstrates the performance of RAR Method IV, which handles elementwise i.i.d. 
bounded CSI errors with unknown distribution. We test the method using elementwise i.i.d. uniform CSI 
errors, where the real and imaginary parts of all are independent and uniformly distributed on [— e, e] 
with e > 0. The probabilistic SOCP method in [21] also has a version for i.i.d. uniform CSI errors and 
is included in our simulation. 

The simulation settings are: Nt = K = 3, p = 0.1, e = 0.02, and the feasible realizations pick-up 
point at 7 = 7dB. The result, presented in Figure 7, illustrates that RAR Method IV provides much better 
performance than the probabilistic SOCP method. Table IV shows the ratios of getting a rank-one RAR 
solution, where we see clearly that encountering high rank RAR solutions is rare. 

VII. Conclusions 

Motivated by the presence of CSI errors in practical systems and the need to avoid substantial SINR out- 
ages among users, we studied a probabilistic SINR constrained formulation of the transmit beamforming 
design problem. Although such formulation can safeguard each user's SINR requirement, it is difficult to 
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Fig. 7. Performance under i.i.d. uniform CSI errors. N t = K = 3; p = 0.1; e = 0.02. 

process computationally due to the SINR outage probability constraints. To circumvent this, we proposed a 
novel relaxation-restriction (RAR) approach, which features the use of semidefmite relaxation techniques, 
as well as analytic tools from probability theory, to produce efficiently computable convex approximations 
of the aforementioned probabilistic formulation. One of our main contributions is the development of three 
methods — namely, sphere bounding, Bernstein-type inequality, and decomposition — for processing the 
probabilistic SINR constraints. Our simulation results indicated that the proposed RAR methods provide 
good approximations to the probabilistic SINR constrained problem, and they significantly improved upon 
existing methods, both in terms of solution quality and computational complexity. 

At the core of our technical development is a set of tools for constructing efficiently computable convex 
restrictions of chance constraints with quadratic uncertainties. An interesting future direction would be to 
apply these new tools to other transmit beamforming formulations, such as those arising from the frontier 
cognitive radio and multicell scenarios, or perhaps even other signal processing applications. 
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